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ABSTRACT 

We present an analysis of a new, large sample of field blue-straggler stars (BSSs) in the thick 
disk and halo system of the Galaxy, based on stellar spectra obtained during the Sloan Digital Sky 
Survey (SDSS) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE). Using 
estimates of stellar atmospheric parameters obtained from application of the SEGUE Stellar Parameter 
Pipeline, we obtain a sample of some 8000 BSSs, which are considered along with a previously selected 
sample of some 4800 blue horizontal-branch (BHB) stars. We derive the ratio of BSSs to BHB 
stars, Fbss/bhbj as a function of Galactocentric distance and distance from the Galactic plane. The 
maximum value found for Fbss/bhb is 4.0 in the thick disk (at 3 kpc < |Z| < 4 kpc), declining to 
on the order of ^ 1.5 — 2.0 in the inner-halo region; this ratio continues to decline to ^ 1.0 in the 
outer-halo region. We associate a minority of field BSSs with a likely extragalactic origin; at least 5% 
of the BSS sample exhibit radial velocities, positions, and distances commensurate with membership 
in the Sagittarius Stream. 

Subject headings: Galaxy: halo—methods: spectroscopy—methods: radial velocities—stars: blue 
stragglers—stars: blue horizontal branch—stars: binaries 


1. INTRODUCTION 

Blue-straggler stars (BSSs) lie brighter and blueward 
of the main-sequence turnoff region in color-magnitude 
diagrams of globular and open clusters. This is a re¬ 
gion where, if the BSSs are single stars, they should have 
already evolved away from the main sequence. Since 
they lie blueward of the cluster turnoff point, and ap¬ 


pear to linger (or straggle behind already evolved stars), 
they have been traditionally referred to as blue stragglers 
(Sandage 1953). BSSs have been found in all stellar pop¬ 
ulations and environments, including: Population I stars 
in the field (Carney et al. 2005), in open clusters of all 
ages (de March! et al. 2006), and among thick-disk stars 
(Fuhrmann & Bernkopf 1999), as well as among Popula¬ 
tion H stars in globular clusters (Piotto et al. 2004), the 
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Galactic halo (Preston et al. 1994), the Galactic bulge 
(Glarkson et al. 2011), and dwarf galaxies in the Local 
Group (Momany et al. 2007). 

Since their discovery, a large number of studies and re¬ 
views have been published discussing the detection of 
BSSs, their properties, and various hypotheses to ac¬ 
count for the blue-straggler phenomenon (Leonard 1989; 
Stryker 1993; Bailyn 1995; Leonard 1996; Preston 2014). 
Two mechanisms for their origin have been most promi¬ 
nently mentioned: (i) Gollisions between stars, resulting 
in a more massive and apparently younger star (Hills & 
Day 1976); and (ii) coalescence and/or mass-transfer pro¬ 
cesses in binary systems, initially proposed by McCrea 
(1964). The latter scenario appears to be the main chan¬ 
nel for the formation of BSSs, where the more-massive 
star in a binary system transfers material, during its post 
main-sequence stage, to its presently-observed compan¬ 
ion (Boffin et al. 2014). The companion star, now more 
massive, appears anomalously younger. 

Fusi Peed et al. (1992) were among the first to propose 
that BSSs in different environments may have different 
origins. In globular clusters, both of the aforementioned 
processes could act simultaneously. Gollisions are ex¬ 
pected to dominate in the central regions, while in the 
outer regions, BSSs mass transfer dominates. As shown 
by Mapelli et al. (2004), the specific frequency of BSSs 
in globular clusters appears highly peaked at cluster cen¬ 
ters, and rapidly decreases at intermediate radii; it rises 
again in the outer regions of most globular clusters. 

A constrasting trend is found when comparing the fre¬ 
quency of BSSs in different environments (Preston et al. 
1994; Momany et al. 2007); denser stellar environments 
(globular clusters) exhibit BSS frequencies that are lower 
than in less-dense stellar environments (open clusters and 
the Galactic field). One explanation for this observation 
suggests that multiple systems (presumably the progen¬ 
itors of BSSs) can be more easily disrupted when sub¬ 
jected to gravitational interactions within the clusters 
(Preston & Sneden 2000). Even the stars found in the 
low-density outer regions of globular clusters may be sub¬ 
ject to stellar encounters if their intra-cluster motions 
carry them into the central high-density regions. 

Preston & Sneden (2000) argued that field BSSs were 
created primarily by mass-transfer processes. They 
noted that field BSSs are much more common per unit 
luminosity than BSSs in the lower-density outer regions 
of globular clusters. Binary systems, even in those low- 
density regions, may have already been destroyed, which 
would reduce the number of BSSs. Preston & Sneden 
(2000) obtained high-resolution spectroscopy for 62 of 
the 84 so-called blue metal-poor (BMP) stars brighter 
than H ^ 15 from Preston et al. (1994). These stars 
exhibit main-sequence gravities and colors bluer than 
most metal-poor globular cluster turnoffs, 0.15 < {B — 
V)q < 0.35. Preston & Sneden (2000) concluded that 
over 60% of their sample are binaries, and that at least 
50% of them are BSSs. These values suggest that the 
BSSs can be a strong indicator of the presence of bi¬ 
nary systems, and can be used as probe of the binarity 
(in particular close binaries) of the stellar populations in 
the Galaxy. 

Glaspey et al. (1994) obtained moderate-resolution 
spectra of the Li 6707 A doublet for a sample of likely 


held BSSs, and obtained clear evidence for Li depletion, 
concluding that this depletion is a general property as¬ 
sociated with blue stragglers. These authors argued that 
the origin of the Li depletion in their sample stars was 
likely related to large-scale mixing, presumably associ¬ 
ated with the formation of blue stragglers. 

Ryan et al. (2001) and Ryan et al. (2002) have also 
argued in favor of the binarity of Held blue stragglers, 
based on consideration of the so-called ultra-Li-dehcient 
halo stars, proposing that these stars may have had their 
surface abundances of lithium reduced by the same mech¬ 
anism that produces Held blue stragglers. These ultra- 
Li-deHcient halo stars were shown to possess more rapid 
rotation than other ancient metal-poor stars, with the 
spin-up likely due to the mass-transfer process, and/or 
chemical peculiarities associated with mass transfer from 
an AGB progenitor. The Ryan et al. sample of stars were 
cooler than the main-sequence turnoff, hence they were 
referred to by these authors as “blue stragglers to be.” 
These stars thus could be viewed as the low-mass coun¬ 
terpart of Held blue stragglers (Ryan et al. 2002). Fur¬ 
thermore, the process of mass transfer that appears to 
be responsible for creating BSSs must act in cooler stars 
as well. Rocha-Pinto et al. (2001) identify cooler blue- 
straggler candidates (yellow or red stragglers). These 
were referred to as chromospherically young and kine¬ 
matically old stars (GYKOS). 

It is worth noting that the BSS frequency in 
low-luminosity dwarf galaxies and open clusters 
(Nbss/Nbhb^~ 2.2; Momany et al. 2007) is different 
than that derived for the Galactic halo (Nbss/Nbhb 
^ 4; Preston et al. 1994). However, as pointed out by 
Momany et al. (2007), the latter value was found based 
on a sample of only 62 BMP stars distributed along 
different lines of sight, located at different distances, 
and for which no observational star-by-star (BSS/BHB) 
correspondence could be deHned. 

Furthermore, Held-star samples have introduced an¬ 
other puzzle, one at least as intriguing as the blue strag¬ 
gler phenomenon itself. Preston et al. (1994), Preston & 
Sneden (2000) and Carney et al. (2001, 2005) argued that 
binary-star evolution and mass transfer is one path for 
the creation of blue stragglers, and apparently the most 
common among Held stars. However, some metal-poor 
main-sequence Held stars that are hotter than globular- 
cluster main-sequence turnoffs are apparently not bina¬ 
ries (see Preston et al. 1994, for further details). These 
authors accounted for their existence as the debris from 
accreted dwarf satellite galaxies whose star formation 
continued over an extended period. 

This paper aims to examine the issues previously raised 
concerning Held BSSs by gathering a large new sample of 
BSS candidates, in order to improve estimates of the spe- 
ciHc frequency of Nsss/A^bhe (Fbss/bhb) in the Galactic 
Held, and to test the hypothesis of the extragalactic ori¬ 
gin for some BSSs. This paper is outlined as follows. 
Section 2 describes the assembly of the database used to 
identify our sample of Held BSSs and BHB stars, and the 
photometric restrictions adopted for both samples. Sec¬ 
tion 3 summarizes our analysis of the selected spectra, 
and the techniques used to partition them into these two 

^ Nbss 3-ud Nbhb ure the number of BSSs and BHB stars found 
in these environments, respectively. 
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groups. Section 4 derives estimates for the absolute mag¬ 
nitudes and distances for both BSSs and BHB stars, and 
evaluates their frequencies in the Galactic field. Section 5 
presents a kinematic analysis of the field BSS sample, in 
an attempt to associate their radial velocities with stel¬ 
lar streams. A brief discussion and our conclusions are 
provided in Section 6. 

2. PHOTOMETRIC SELECTION AND STELLAR 
PARAMETERS 

2.1. SDSS - DR8 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) 
is a large photometric and spectroscopic survey that cov¬ 
ers about one-quarter of the northern sky. This paper 
makes use of the eighth data release (DR8), which con¬ 
tains the data obtained by the survey through January 
2011 (Aihara et al. 2011). The first restrictions adopted 
in the selection of field BSSs and BHB stars follow the 
same criteria used by Xue et al. (2008). 

2.2. Color Cuts 

Blue-straggler and blue horizontal-branch candidates 
are often identified in the SDSS color-color diagram re¬ 
gion with the following color cuts: 0.60 < [u — g) < 1.60 
and —0.50 < {g—r) < 0.05 (Xue et al. 2008; Deason et al. 
2011). This region is preferentially populated by A-type 
stars. For this work, we first obtained a sample of 19771 
objects within these color limits, and with available spec¬ 
tra, but with no restrictions on apparent magnitudes or 
location on the sky within the SDSS footprint, and cov¬ 
ering the range 14.5 < g < 21, the bright limit being 
set by saturation in the SDSS imaging scans. Note that 
the ugriz magnitudes and color indices are all corrected 
for absorption and reddening using the extinction maps 
from Schlegel et al. (1998). 

Previous efforts to study stars in this color region usu¬ 
ally targeted only BHB stars (Sirko et al. 2004; Xue et al. 
2008; Deason et al. 2011; Xue et al. 2011). In this work, 
we are interested in selecting both BSS and BHB candi¬ 
dates. However, photometry alone does not uniquely dis¬ 
tinguish main-sequence A-type stars (of which the metal- 
poor stars are BSS candidates) from giant (BHB) can¬ 
didates. Furthermore, in this color region, the A-type 
stars can be mixed with stars of cooler spectral types, 
such as metal-poor F- and G-type stars. Thus, we also 
made use of available spectroscopy to trim our sample, 
as described below. 

2.3. Stellar Parameters 

The spectra of the selected sample have stellar param¬ 
eters determined by the current version of the SEGUE 
Stellar Parameter Pipeline (SSPP, Lee et al. 2008a,b; 
Allende Prieto et al. 2008; Smolinski et al. 2011; Lee 
et al. 2011). Erom the SSPP, we obtained values of effec¬ 
tive temperature (Teff), surface gravity (logg), metal- 
licity ([Ee/H]^), and heliocentric radial velocity (Uhrv)- 
To select the best BSS candidates (as opposed to fore¬ 
ground A-type stars in the disk populations), we applied 

^ The abundance of an element in relation to another is defined 
in the usual manner, using the notation: [A/B] = log(Nj\^/NB)* — 
log[N/^^/N b)q^ where and Nq are the number of atoms of 
elements A and B in the star (★) and in the Sun (O), respectively. 



Fig. 1.— Normalized spectrum of a BHB star (red line) and 
a BSS (black line) for the same T^ff (~8100 K) in the 
region. The dashed line is drawn at 20% below the continuum, 
highlighting the difference between the Balmer-line wings. 

a restriction in metallicity, [Fe/H] < —0.4, reducing the 
sample to 18560 objects. 

It is relatively straightforward to identify BSSs in glob¬ 
ular clusters because the horizontal-branch stars, whose 
colors blend with hotter main-sequence stars, are seen at 
different apparent magnitudes (BHB stars are typically 

2 magnitudes brighter than BSSs, as noted by Preston 
et al. 1994; Deason et al. 2012). Unfortunately, these 
two populations cannot be easily separated in the field, 
because the stars are not at the same distance. 

In the temperature range of A-type stars, the Balmer 
lines are the main tool to separate A-type main-sequence 
from giant stars (BSSs from BHB stars, respectively). 
The temperature modifies the depth of the Balmer lines, 
while the surface gravity changes the width of its wings 
(as initially noted by Pier 1983). These effects have 
a clear impact on the appearance of the spectrum, as 
shown in Figure 1: for stars with the same Teff (in this 
case, ~8I00 K), the Balmer-line wings are broader for 
higher surface gravities. 

3. THE BLUE-STRAGGLER AND BLUE 
HORIZONTAL-BRANCH SAMPLES 

In this section we describe the methods used to sepa¬ 
rate BSSs from BHB stars, using the available medium- 
resolution SDSS/SEGUE spectroscopy. Deason et al. 
(2012) states that the differences between BSSs and BHB 
stars becomes most evident for spectra with S/N>5 (see 
the bottom panels in Figure 3 of Deason et al. 2012). For 
this work, we prefer to only use stars with spectra having 
an average signal-noise ratio (over the optical spectrum) 
of (S/N) > 10. This limit allows for a confident distinc¬ 
tion between BSSs and BHB stars where they are photo¬ 
metrically indistinguishable, and for the reliable stellar 
parameters from the SSPP. 

Using the parameters calculated by the SSPP, we can 
compare the surface gravities and effective temperatures 
of BSSs with BHB stars (similar to Wilhelm et al. 1999; 
Deason et al. 2012). However, the atmospheric parame¬ 
ters calculated for stars with temperatures above 7500 K 
(where the BSSs and BHB stars are usually found) re¬ 
quire additional quantitative validation. This is accom- 
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plished by comparing the results of the SSPP with other 
methods for separating the BSSs from BHB stars, which 
generally required higher S/N spectra, typically greater 
than (values typically achieved for SDSS spectra of 
stars with apparent magnitudes g < 18). 

It is also important to note that, even though stars 
with g > 18 are usually not suitable for analysis by 

all the methods presented below, they provide valuable 
information on the populations of BSSs and BHB stars 
(especially since, as shown below, the majority of likely 
BSSs in our sample are fainter than g = 17). For the 
following analysis, we divided the sample in two: (i) 
The Bright Sample {g < 18) and (ii) the Faint Sample 
(5 > 18). 

3.1. The Bright Sample: 14.5 < g < 18 

Since the Bright Sample contains stars having spectra 
with higher S/N ratios, it is possible to apply a series of 
quantitative tests to distinguish between BSSs and BHB 
stars. These are described below. 

3.1.1. SSPP Restrictions 

The parameters provided by the SSPP were used to 
identify BSSs and BHB stars. Since the Balmer lines are 
strong in the 7500 K < < 10000 K range, it is pos¬ 

sible to effectively distinguish between BSSs and BHB 
stars from their loggdistribution (Wilhelm et al. 1999). 
Figure 2 summarizes the adopted restrictions. The line 
separating these two groups of A-type stars is log g = 
3.80, as estimated from the logg distributions, for stars 
with g < 18 (see top panel of Figure 5). This value 
represents a threshold at which, statistically, the con¬ 
tamination of the BSS population by BHB stars is mini¬ 
mal. The limits lie at Scr from the peaks in the top panel 
of Figure 5 centered at log 5 = 3.38 (cr = 0.14 dex) and 
logg= 4.33 {a = 0.18 dex), which represent the BHB 
stars and BSSs, respectively. 

3.1.2. The fm vs. Do .2 Method 

The parameters used in the following analysis are ob¬ 
tained through fits to a Sersic profile (Sersic 1968), 
which describes the shape of the Balmer spectral features 
(Clewley et al. 2002): 

S'(a;) = n - a • exp[-(|A - Ao|/6)“], (1) 

where, n, a, 6, and c are free parameters. The param¬ 
eter n fits the continuum level (~1 for our normalized 
spectra), Aq is the wavelength at the line center, and 
a defines the scale-depth parameter. The parameters b 
(scale width) and c (scale shape) will be subsequently ex¬ 
amined. They are used in the scale width-shape method 
(Clewley et al. 2002), as described in Section 3.1.3. 

There is no fundamental relation between a and the 
other Sersic parameters, but the a values have well- 
defined Gaussian distributions around a mean value for 
each Balmer line. For our sample, these average val¬ 
ues are: (a^) = 0.709, (a.y) = 0.694, and (a^) = 0.769. 
We refined the estimated Sersic profile parameters by re¬ 
adjusting the function with (a) fixed for each line. 

The fm versus D 0.2 method was first proposed by Pier 
(1983); it compares Teff and log g indirectly, by measuring 
specific quantities in the stellar spectra. The parameter 
fm is the deepest flux relative to the continuum level, 



9000 8000 7000 6000 5000 


Teff(K) 

Fig. 2. — Distribution of the SSPP Teff vs. loggfor stars with g 
< 18. There are two stellar groups concentrated in Teff > 7500 K: 
one for log^'^ 3.4 (BHB stars: red dots) and other with logg~ 
4.3 (BSSs: black dots). The gray squares are cooler stars and 
are not considered further in our analysis. For convenience of 
the reader, we highlight the positions where the RR Lyrae stars 
are often found using the blue rectangle (Wilhelm et al. 1999), 
where the gap with no stars present is the expected location of the 
instability strip; the purple rectangle shows the region where the 
occupied by main-sequence turnoff stars (MSTO). Overplotted is 
a Yale-Yonsei isochrone (Demarque et al. 2004) for 12 Gyr and 
[Fe/H]= -2.0. 

usually measured at the center of the line (it is sensitive 
to Teff ). The parameter D 0.2 is the width of the line mea¬ 
sured at 20% below the adjusted continuum level (Beers 
et al. 1992; Sirko et al. 2004), and provides an indicator 
of logy. As seen in Figure 1, when logy is higher, the 
greater D 0.2 becomes. To obtain D 0 . 2 , we subtracted the 
roots of the inverse function when S{x) = 0.8n, following 
Clewley et al. (2002), which yields: 

no ,2 = 2 • 6 • [log(n/(5 • (a))]b (2) 

We adopted the H^ feature to compare the parameters 
fm and Do .2 for this method, indicated by 5 in the sub¬ 
script of the derived parameters (see Sirko et al. 2004, for 
further details). Figure 3 shows the behavior of Do, 2 S, as 
a function of fmS, for stars with g < 18. The stars con¬ 
centrated at {fmS: Do. 2 s) (0.26; 25) are the BHB stars; 
the stars with larger Dq ,25 are BSSs. The remaining stars 
have cooler spectral types. 

To compare the restrictions adopted by the SSPP with 
the limits in the fm vs. D 0.2 method, we identified the 
stars in Figure 3 according to the regions described in 
Figure 2; black dots represent BSSs, red dots are BHB 
stars, and the gray squares are cooler stars. Note that 
these divisions exhibit good agreement with the separa¬ 
tion based on the SSPP stellar-parameter values (~92% 
of stars satisfy both restrictions). 

3.1.3. The Scale Width-Shape Method: c vs. b 

An additional method often adopted to separate BSSs 
from BHB stars is the scale width-shape method (c vs. 
b), first described by Clewley et al. (2002). This method 
separates these stars through the parameters set by the 
Sersic profile (with the parameter a fixed). In order to 
apply this approach, we calculated the average of the 
parameters c and b for the H^, H.^ and H^ lines, fol- 
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Fig. 3.— Comparison of vs. Do.25? color-coded according to 
the restrictions adopted by the SSPP (Figure 2). The distributions 
of BSSs and BHB stars exhibit some overlap beyond the adopted 
limits (dashed lines). The adopted criteria to select BSSs are: (i) 
0.15 < W < 0.27, Do.25 > 29; (ii) 0.27 < < 0.40, Do.25 > 

—90 • fni 5 + 53. Accordingly, the criteria for BHB candidates are: 
(i) 0.15 < fms < 0.27 and 17 < Dq 25 < 29; (ii) 0.27 < W < 0.40 
and 17 < Do .25 < -90 ■ W + 53. 

lowing Deason et al. (2011). The average of these pa¬ 
rameters is shown in Figure 4, where and are 
color-coded by the adopted regions described in Figure 2. 
As can be seen, there is good agreement when comparing 
this method with separation based on the SSPP stellar- 
parameter values (~95% of stars satisfy both restric¬ 
tions). We note that the parameters b and c could be 
used, in principle, by themselves to separate BSS from 
BHB (see Preston 2014). The mixed stars that passed 
the initial color cuts now exhibit a clear separation in 
the Ci 3 js versus bp-ys plane; the gap between them dis¬ 
tinguishes between BSSs and BHB stars, as defined by a 
fourth-order polynomial: 

bp^s = -1-9.20 — 46.32 • (cpys) + 82.24 • {cp^s)^ 

-23.36 ■ {cp^sf - 10.82 ■ (cp^s)^. (3) 

To select BSSs with g < 18, we combined the SSPP 
restrictions, /m vs. II 0 . 2 , and the scale width-shape 
method. Whereas the constraints of the SSPP proved 
fully compatible with the methods based on the analysis 
of the Balmer lines, we further restrict the BSS and BHB 
candidates according to the limits: BSSs have log g >3.8 
and the BHB stars lie between 3.0 < \ogg < 3.8. The 
mean error on effective temperature given by the SSPP 
is about 150 K (Lee et al. 2008b), which is slightly higher 
than the average error of the parameters for our sample 
(~130 K). We selected only the stars with > 7650 K. 
The stars satisfying these cuts are within 3cr of the mean 
of both peaks in logg , as shown in the top panel of Fig¬ 
ure 5, resulting in 4838 BSSs and 4380 BHB stars. 

3.2. The Faint Sample: 18 < g < 21 

The traditional methods (/m vs. D 0 . 2 , and c vs. b) 
cannot properly evaluate the stars in our sample with 
g > 18, due to the lower (S/N) ratios in SDSS/SEGUE 
spectra of these fainter stars (Clewley et al. 2002; Sirko 
et al. 2004). Thus, for the 5663 stars in the Faint Sample, 



Fig. 4. — Comparison between cpjg and bp-ys, with color-coding 
according to the restrictions adopted by the SSPP. As verified 
with the /m vs. D 0.2 method, the limits adopted by the SSPP 
for separating the BSSs and BHB stars agree very well (~95% of 
stars satisfy both restrictions). The limit line that separates these 
two groups is defined by the fourth-order polynomial given by 
Equation 3. 

we employed only the SSPP parameters to separate the 
BSSs and BHB stars. 

The log g distribution for the Faint Sample is shown in 
the lower panel of Figure 5. The peak centered at log g = 
3.57, which represents BHB stars, has a high dispersion 
(cr = 0.33 dex). This increases the size of the overlap 
region between the two distributions and, if we were to 
consider the Scr region as a limit, it would exclude almost 
all BSSs, which are peaked at log 5 = 4.27 (tr = 0.20 dex). 

As the stars analyzed in the lower panel of Figure 5 
are also well-described by Gaussian distributions, we as¬ 
sessed the probability of finding BHBs using the peak 
positions and values up to 3 (t away from them. The 
adopted restriction for the BHB distribution is log g < 
3.66 (3cr away from the BSSs peak). 

Since our goal is to select BSSs candidates, we adopt a 
more relaxed criteria to its distribution, even though it 
leads to greater contamination from the BHB population. 
Estimating the amount of stars coming from each distri¬ 
bution, we analyzed the contamination for every sigma. 
For example, there are 3646 stars with log(;> 4.27 (lo¬ 
cated at log g BHB + 2(Tbhb), but the probability of find¬ 
ing BHB stars in this range of log g is approximately 2%. 
Thus, stars with log( 7 > 4.27 have a 98% probability of 
being BSSs. 

We conclude that, for g > 18, the most reliable BHB 
candidates are found in the surface-gravity range 3.00 
< log 5 < 3.66 (416 objects). Accordingly, for BSSs we 
adopted the limit of log g > 3.92, which is Icr away from 
the BHB peak, with an average probability of ~ 88 %. 

Gombining the Bright and Faint Samples, we obtained 
a total of 8001 BSS candidates, with an average proba¬ 
bility of ~95%. We also obtained 4796 BHB candidates, 
with an average probability of correct selection of ~98%. 

Figure 6 shows the distributions of log g from the SSPP 
for our full sample of BSSs and BHB stars, in intervals of 
apparent magnitude. As is clear from inspection of this 
sample, the higher gravity BSSs begin to dominate over 
the lower gravity BHB stars for g > 17. Figure 7 shows 
the distribution of derived [Fe/H] for these same stars. 
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Fig. 5.— Top panel: Distribution of log ^ for stars with g < 18. 
The peak centered at \ogg= 3.38 (cr = 0.14 dex) represents the 
distribution of BHB stars; that at \ogg= 4.33 (cr = 0.18 dex) 
applies to BSSs. Bottom panel: Distribution of log ^ for stars 
with g > 18. The peak centered at log 51 = 3.57 (cr = 0.33 dex) 
represents the distribution of BHB stars; that at log 5 '= 4.27 {cr 
= 0.20 dex) represents BSSs. It can be seen that the fraction of 
BHB stars decays quickly at high surface gravity for the fainter 
stars. 
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Fig. 7. — [Fe/H] distributions for our sample stars (left panels: 
BSSs; right panels: BHB stars) for the Bright Sample (upper 
panels) and the Faint Sample (bottom panels). The gray lines 
indicate the mean of each [Fe/H] distribution. 


SDSS ugriz system (e.g., Kinman et al. 1994; Preston 
et al. 1994; Preston & Sneden 2000; Carney et al. 2001, 
2005). In this work, we used the transformations of Zhao 
& Newberg (2006), which were derived from SDSS stars 
with known UBVRI photometry, including a number of 
BSSs and BHB stars: 


Vo = g - 0.561 -{g-r)- 0.004, (4) 


and 
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Fig. 6. — Histograms of logg, evaluated in different intervals of 
g-band apparent magnitude. The number of BSSs dominates over 
the number of BHB stars at fainter apparent magnitude (g > 17). 


in this case for the Bright Samples and Faint Samples of 
BSSs and BHB stars. Although the mean metallicities of 
the BSSs and BHB stars are similar, the BSSs are more 
prevalent at higher metallicity than the BHB stars. 

4. DISTANCE CALIBRATIONS AND FREQUENCIES 

In this section we obtain an absolute magnitude cali¬ 
bration for our full sample of 8001 BSSs and 4796 BHB 
stars, extending over the magnitude range 14.5 < g < 21. 

A number of studies adopt calibrations for apparent 
magnitudes based on photometric systems other than the 


(H-V)o = 0.916-( 5 -r-)-b 0.187. (5) 

According to Beers et al. (2012), the valid range for 
these conversions is —0.5 < g — r < 1.0, consistent with 
the color indices adopted in this work. In addition, the 
calculated values are within the limits proposed by Car¬ 
ney et al. (2005) to define field BSSs. 

4.1. Absolute Magnitudes of Field BSSs 

Based on the V-band magnitude, it is possible to esti¬ 
mate the absolute magnitude for BSSs using the relation 
from Kinman et al. (1994): 


Mvbss = 1-32 + 4.05 • {B - V)o - 0.45 • [Fe/Hj. (6) 

This relation was first determined from BSSs in glob¬ 
ular clusters by the work of Sarajedini (1993, 1994). As¬ 
suming {B — V)o = 0.1 and [Fe/H] = —1.7 as typical 
values for metal-poor A-type stars selected in this work, 
we verify that the field BSS candidates have My ~ 2.5. 
This value is in agreement with blue stragglers observed 
in globular clusters and dwarf galaxies (Sarajedini 1993; 
Momany et al. 2007). 

Deason et al. (2011) also developed a calibration for the 
g-band absolute magnitude for field BSSs. They made 
use of stars in Stripe 82 that are members of the Sagit¬ 
tarius Stream, combining the apparent magnitudes and 
the results obtained by Watkins et al. (2009). The latter 
used RR Lyrae stars to estimate the distance to stars in 
the stream in the right ascension range 25° < a < 40°, 
with Dsgr = 26.1 ± 5.6 kpc. The absolute-magnitude 
calibration in the g-band from Watkins et al. (2009), 
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MgBss = 3.108 + 5.495 • {g — r), is similar to that pre¬ 
sented by Kinman et al. (1994) for [Fe/H] = —1.5, and 
also agrees with the mean metallicity of stars identified 
in the Sagittarius Stream ([Fe/H] = —1.43, Watkins 
et al. 2009). However, unlike the relation used by Kin- 
man et al. (1994), Deason et al. (2011) found no strong 
dependence on metallicity. This is probably due to the 
use of BSSs in the Sagittarius Stream, which already have 
a mean [Fe/H] = —1.4, and also have similar luminosities 
(Watkins et al. 2009). The adopted error for this cali¬ 
bration is (iMg = 0.5 mag, primarily arising from the 
expected error in the distance to the Sagittarius Stream. 

The distances are systematically smaller for the 
method calibrated to BSSs in the Sagittarius Stream, but 
both relations agree within 2cr. We adopt the calibration 
given in Equation 6 of Kinman et al. (1994), because it 
was determined from BSSs in globular clusters with dif¬ 
ferent metallicities covering a wide range of colors and 
absolute magnitudes. Furthermore, we take into account 
that the error in the absolute magnitude of each star is 
certainly greater for the fainter objects. 

4.2. Absolute Magnitudes of Field BHBs 

We employed two different calibrations to determine 
absolute magnitudes for the BHB sample: (1) Deason 
et al. (2011) derived Mg using SDSS photometry of 10 
globular clusters from An et al. (2008); (2) Fermani 
& Schonrich (2013) proposed a calibration with a de¬ 
pendence on the metallicity, arguing that the relation 
adopted by Deason et al. (2011) is systematically fainter 
for metal-poor stars and brighter for metal-rich stars. 
The absolute magnitudes for Deason et al. 2011 (Mgo) 
and Fermani & Schonrich 2013 (Mgps) are given by: 

MgD = 0.434 - 0.169 ■ [g - r) + 2.319 • {g - rf 

-f20.449 • (g - rf + 94.517 ■ {g - rf, (7) 

MgFs = 0.0075 • e-i4 °49-D 
-H0.04-([Fe/H]-H 3.5)2-h 0.25. 

The range of Mg for BHB stars in globular clusters 
evaluated by Deason et al. (2011) is 0.45 < Mg < 0.65, 
and their metallicities are in the range [—2.3 : —1.3]. 
The authors did not find any significant dependence on 
metallicity for Mg. 

According to Fermani & Schonrich (2013), the lack of 
metallicity dependence on the absolute magnitude cali¬ 
bration of Deason et al. (2011) is primarily due to un¬ 
certainties in their adopted distance moduli, reddening, 
and presumed systematic differences between clusters 
and field halo stars; they thus argued that Deason et 
al. could not rule out a detectable metallicity depen¬ 
dence on field BHB stars. For the distance calculations 
presented below, we explore both Mg calibrations. 

4.3. The Relative Frequencies of BSSs and BHB Stars 
4.3.1. Galactocentric Distances 

In order to determine the ratio between BSSs and BHB 
stars (Fbss/bhb); we first proceed by counting them di¬ 
rectly, as a function of their distance from the Galac- 



Rdi, (Kpc) 



Roll (kpc) 



F*FS13 “ *^D11 


Fig. 8.— (a): Comparison between distances estimated using 
the absolute magnitude calibrations of Deason et al. (2011) and 
Fermani & Schonrich (2013). The solid line represents itpsis = 
Run, the red dashed line is a linear fit to the data, and the gray 
region is the ± 3 (t limit of the data. Panels (b) and (c) show, 
respectively, the residuals for the distribution and its dispersion, 
represented by a Gaussian fit. The dark gray region represents 
the limit within Icr of the average value (~0.34 kpc) tor the 
distribution; the light gray region is the limit within 2a of the 
average value. 

tic center. The Galactocentric distance, R, is calculated 
based on the distance, D, from the Sun (in kpc): 

= io[(™-m)/5-2]^ (9) 

and 

= (_Rg — D cosbcosl)^ -|- (D sin 6)2 -|- (T)cos6sin/)2, 

where 6 and I are the Galactic latitude and longitude, 
respectively, and Rq is the distance of the Sun to the 
Galactic center. We adopted Rq = 8.5 kpc to be con¬ 
sistent with the recent kinematic analysis from Garollo 
et al. (2010). 

Figure 8 compares the Galactocentric distances ob¬ 
tained using the BHB calibrations from Deason et al. 
(2011) (i?Dii) and Fermani & Schonrich (2013) (i^FSis)- 
The i?FSi 3 distances are, on average, ^ 2.5% larger than 
7 ?dii, which is within the errors for both calibrations. 
Also, from panels b and c of Figure 8, it is possible to see 
that, even though there is a larger spread for distances 
above 40 kpc, there are no trends between the calibra¬ 
tions. The zero-point offset for this comparison is 0.27 
kpc, with a standard deviation of 0.39 kpc. 

It is clear that the differences between these two ap- 
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Fig. 9.— Comparison of the “as observed” Fggs/BjjBi ^ 
function of Galactocentric distance, using the absolute-magnitude 
calibrations from Deason et al. (2011) and Fermani & Schonrich 
(2013). The gray line represents the value of Fggg/gjjg estimated 
in the Solar Neighbourhood by Preston et al. (1994), and the 
dashed line represents the value of Fggg/gjjB found in nearby 
dwarf galaxies by Momany et al. (2007). Corrections for the 
sampling biases of these populations are described in the text. 

proaches increase for larger distances. For this reason, 
we compared the “as observed” frequency FBsg/BHB for 
both Rdii and Rfsi 3 j to gauge the possible impact of 
choosing one or the other of these calibrations. Results 
are shown in Figure 9 (values can be found in Table 1). 
Both frequencies agree within 3cr, and the residual val¬ 
ues have an average of zero, with a standard deviation 
of 0.16 kpc. Thus, we choose to employ the calibration 
from Deason et al. (2011) to represent the distances to 
BHB stars in the Galactic field, since we found that the 
introduction of the [Fe/H] as a parameter does not have 
any impact on the derived distances. 

For completeness, Figure 10 shows the derived dis¬ 
tances of the BSSs and BHB stars in our sample plotted 
in the XYZ planes of the Galaxy. As expected, the BHB 
stars explore to greater distances than the BSSs (and are 
undersampled in the Solar Neighborhood relative to the 
BSSs, due to the SDSS bright limit), so care must be 
taken to account for these limitations in obtaining our 
Hnal derived relative frequency estimates. 


Fig. fO.— Distribution of our sample stars (left panels: BSSs; 
right panels: BHB stars) in the XZ plane (upper panels), and in 
the XY plane (the Galactic plane; lower panels). The gray lines 
are drawn for reference to include the location of the Galactic 
center. 



R (kpc) 

Fig. 11.— Distribution of height above the Galactic plane, Z, 
as a function of Galactocentric distance, R, for the BSSs (black 
circles) and BHB stars (red squares) in our sample. The blue 
rectangles limit the regions where we have obtained the relative 
stellar counts, avoiding the difficulties (under/over counts of 
BSSs and BHB stars at the bright/faint limits of our sample), 
as noted in the text: 3 kpc < |Z| <30 kpc and 6 kpc < R < 30 kpc. 


4.3.2. Refined Relative Frequency Estimates 

The frequencies presented in Figure 9 were derived 
to compare the differences resulting from the adoption 
of different absolute magnitude calibrations. However, 
as foreshadowed above, we cannot evaluate the value of 
Fbss/bhb hi the held without taking into account that 
there exist regions of the sampled volume that are not 
simultaneously occupied by BSSs and BHB stars. Near 
the bright magnitude limit of the SDSS scans {g ^ 14.5) 
the closest selected A-type stars will preferably be BSSs, 
as BHB stars at similar distances will generally be satu¬ 
rated. Indeed, at |Z| < 1 kpc, our sample includes more 
than 250 BSSs and only 4 BHB stars. To be safe, we 
have obtained our counts avoiding the region within |Z| 
= 3 kpc. At larger distances (|Z| > 30 kpc and R > 
30 kpc), our sample preferentially includes BHB stars, 


and only a few BSSs. Figure 11 summarizes the restric¬ 
tions we have adopted to derive the relative frequencies, 
Fbss/bhb, shown in Figure 12. 

From inspection of Figure 12, the relative blue- 
straggler frequency in the Galactic halo system is consis¬ 
tent with the value of F bss/bhb obtained for the nearby 
dwarf galaxies, but not with that obtained for the Solar 
Neighborhood. The average value of Fbss/bhb from R 
= 7 kpc to R = 27 kpc is 1.83 ± 0.23. 

Preston et al. (1994) evaluated the Fbss/bhb in the 
Solar Neighborhood using stars within 2 kpc of the Sun, 
a region our data precludes us from exploring. However, 
we can also examine the Fbss/bhb values from another 
perspective, evaluating the change in frequency with dis¬ 
tance from the Galactic plane, |Z|, from |Z| = 3 kpc to 
|Z| =30 kpc, still considering R between 6 and 30 kpc, as 
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TABLE 1 

Stellar Counts in Galactocentric Distance Intervals 


R(kpc) 

Nbss 

Deason et al. (2011) 

Fermani & Schonrich (2013) 

Nbhb 

Fbss/bhb 


Nbhb 

Fbss/bhb 

^-^BSS/BHB 

06-07 

155 

99 

1.57 

±0.14 

93 

1.67 

±0.14 

07-08 

242 

123 

1.97 

±0.19 

129 

1.88 

±0.21 

08-09 

293 

159 

1.84 

±0.25 

147 

1.99 

±0.24 

09-10 

366 

140 

2.61 

±0.23 

143 

2.56 

±0.26 

10-11 

426 

146 

2.92 

±0.36 

130 

3.28 

±0.35 

11-12 

501 

159 

3.15 

±0.42 

164 

3.05 

±0.49 

12-13 

508 

208 

2.44 

±0.46 

181 

2.81 

±0.44 

13-14 

470 

168 

2.80 

±0.33 

181 

2.60 

±0.40 

14-15 

483 

184 

2.63 

±0.39 

191 

2.53 

±0.36 

15-16 

463 

220 

2.10 

±0.36 

179 

2.59 

±0.35 

16-17 

437 

195 

2.24 

±0.27 

202 

2.16 

±0.36 

17-18 

465 

205 

2.27 

±0.30 

206 

2.26 

±0.28 

18-19 

354 

201 

1.76 

±0.30 

195 

1.82 

±0.30 

19-20 

386 

172 

2.24 

±0.22 

180 

2.14 

±0.23 

20-21 

318 

199 

1.60 

±0.30 

171 

1.86 

±0.28 

21-22 

335 

147 

2.28 

±0.20 

166 

2.02 

±0.23 

22-23 

271 

167 

1.62 

±0.30 

162 

1.67 

±0.26 

23-24 

207 

145 

1.43 

±0.20 

150 

1.38 

±0.21 

24-25 

194 

135 

1.44 

±0.17 

144 

1.35 

±0.17 

25-26 

192 

115 

1.67 

±0.17 

127 

1.51 

±0.16 

26-27 

140 

115 

1.22 

±0.21 

98 

1.43 

±0.18 

27-29 

229 

173 

1.32 

±0.15 

166 

1.38 

±0.17 

29-32 

225 

174 

1.29 

±0.16 

216 

1.04 

±0.17 

32-36 

166 

235 

0.71 

±0.16 

215 

0.77 

±0.13 

36-45 

89 

407 

0.22 

±0.09 

413 

0.22 

±0.10 

Total 

7915 

4391 

1.80 

±0.22 

4349 

1.82 

±0.23 


Note. — There are 76 BSSs and ~65 BHB stars within R < 6 kpc. There are also 10 BSSs and 
~360 BHB stars that lie between 45 kpc < R < 100 kpc. 



Fig. 12. — Comparison of the reative frequency, Fggg/BHB! 
a function of the Galactocentric distance, for all sample (gray 
circles) and considering the limited regions defined in Figure 11 
(black circles). The gray line represents the value of Fegg/gHB 
estimated in the Solar Neighborhood by Preston et al. (1994) and 
the dashed line represents the value of Fbss/bhb found in nearby 
dwarf galaxies by Momany et al. (2007). The Fggg/BHB found 
for the Galactic halo system are more consistent with the value of 
Fbss/bhb nearby dwarf galaxies than to that found for the 
Solar Neighborhood. 

shown in Figure 13. From inspection of this Figure, the 
relative frequency of BSSs increases at distances closer 
to the Galactic plane, reaching a maximum F^sg/gHB 
=4.11 ± 0.47 at 3 kpc < |Z| <4 kpc, in agreement with 
the previous calculation of the Solar Neighborhood fre¬ 
quency. 



Fig. 13.— Comparison of Fggg/BHB? a function of the 
distance from the Galactic plane, |Z|, considering the limited 
regions defined in Figure 11. The gray line represents the value 
of Fggg/BHB estimated in the Solar Neighborhood by Preston 
et al. (1994) and the dashed line represents the value of Fegg/g^B 
found in nearby dwarf galaxies by Momany et al. (2007). The 
relative frequency of BSSs rises closer to the Galactic plane, in 
agreement with the Solar Neighborhood frequency obtained by 
Preston et al. (1994). 


5. POSSIBLE IMPACT OF THE PRESENCE OF STREAMS 

It has been shown by a number of recent investiga- 
tions (e.g., Majewski et al. 2003; Belokurov et al. 2006; 
Koposov et al. 2012) that the relative importance of de¬ 
bris streams increases with distance from the Galactic 
center, hence we wish to explore the possible impact of 
this on our results, making use of our observed stellar 
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radial velocities. 

5.1. Distribution of Galactocentric Radial Velocities 

Due to their relatively low spatial density, BSSs having 
R > 25 kpc cannot be readily distinguished as members 
of stellar streams in the outer region of the Galactic halo 
based on their distances alone. Analysis of their radial 
velocities can, however, be used to check for the possi¬ 
ble presence of “extragalactic” blue stragglers (i.e., those 
captured from dissolved satellites) among the stars in the 
field population. 

In order to carry out this analysis, the heliocentric 
radial velocities (Vhrv) of BSSs were converted to the 
Galactic standard of rest (GSR) frame. The adopted ve¬ 
locity for the Local Standard of Rest (Vlsr) is 220 km 
s“^, and a solar motion {U,V,W) = (-1-10.1,-1-4.0,-1-6.7) 
km s“^, defined in a right-handed Galactic system with U 
pointing towards the Galactic center, V in the direction 
of Galactic rotation, and W towards the north Galactic 
pole (Hogg et al. 2005). Hereafter, Vqsr is the radial ve¬ 
locity in the GSR frame (i.e., the radial velocity compo¬ 
nent along the star-Sun direction, corrected for Galactic 
rotation). Then, 


Vgsr = Vhrv + (Vlsr + 5-2) sin Z cos 5 

-1-10 cos Z cos 5 -I- 7.2 sin 6. (11) 

Using the Vqsr values, one can search for stellar 
streams associated with the various structural compo¬ 
nents of the Galaxy, including the thick disk, the inner 
halo, and the outer halo. By consideration of the differ¬ 
ent distributions, segregated by distance from the Galac¬ 
tic plane, |Z|, the inner-halo component of the Milky 
Way is expected to dominate the population of halo stars 
found at distances up to |Z| ^10-15 kpc. The outer-halo 
component is expected to dominate in regions beyond 
|Z| ^15-20 kpc (Carollo et al. 2007, 2010; Beers et al. 
2012). It is important to note that the limits used to sep¬ 
arate these in the previous work considered Z^ax and not 
|Z|. The derived parameter Z„iax depends on the adopted 
gravitational potential. In this work we assume that, if 
there are BSSs whose origin is extragalactic, differences 
in their radial-velocity distributions in the Galactic stan¬ 
dard of rest will be noticed in intervals of either Z^ax or 
|Z|. 

The velocity distributions of BSSs and BHB stars, con¬ 
sidering the approximate limits where each Galactic com¬ 
ponent becomes dominant, are shown in Figures 14 and 
15, respectively. The mean value for the distributions is 
usually close to zero (with dispersion of about 100-150 
km s“^), representing the typical values of radial veloc¬ 
ities for halo stars (Hogg et al. 2005). The mean value 
of velocities for the BSSs and BHB stars for the most 
distant stars is shifted toward negative values. This con¬ 
tribution is likely due to a stellar stream, such as Sagit¬ 
tarius, in both the distance from the Galactic center (R - 
left panels) and in distance from the Galactic plane (|Z| - 
right panels). The dashed lines represent Gaussian fits of 
these distributions, while the vertical lines indicate their 
mean values, which are always negative for the most dis¬ 
tant stars. These behaviors are evident from inspection 
of Figures 16 and 17. 
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Fig. 14.— Distribution of Vqsr for BSSs, considering the 
regions limited by R (left panels) and |Z| (right panels). The main 
distribution found in each panel (usually with mean value close to 
zero and dispersion of about 100-150 km s~^) represents typical 
values of radial velocities for halo stars. The vertical lines show 
the mean value of Vqsr for each single Gaussian distribution 
fitted (dashed black line). It is evident that there is a shift in the 
mean value of the Vqsr when we reach the outer-halo region seen 
in the bottom panels. 
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Fig. 15.— Distribution of Vgsr for BHBs, considering the 
regions limited by R (left panels) and |Z| (right panels). The main 
distribution found in each panel (usually with mean value close to 
zero and dispersion of about 100-150 km s~^) represents typical 
values of radial velocities for halo stars. The vertical lines show 
the mean value of Vqsr for each single Gaussian distribution 
fitted (dashed black line). It is evident that there is a shift in the 
mean value of the Vqsr when we reach the outer-halo region seen 
in the bottom panels. 


5.2. Comparison with the Sagittarius Stream 

We now compare the Vqsr foi" BSSs and BHB stars 
found at |Z| > 17 kpc with typical values for the Sagittar¬ 
ius Stream (Majewski et al. 2003), in order to verify the 
hypothesis that at least a portion of the BSSs and BHB 
stars in our sample can be associated with this stream, 
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Fig. 16. — Mean values of Vqsr shown in Figures 14 and 15, as 
a function of R. The black circles are the BSSs and the red are 
squares the BHBs. Error bars were estimated using the standard 
deviations of the fitted Gaussians. 
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Fig. 17.— Mean values of Vqsr shown in Figures 14 and 15, 
as a function of |Z|. The black circles are the BSSs and the 
red squares are the BHBs. Error bars were estimated using the 
standard deviations of the fitted Gaussians. 




Fig. 18.— Galactic celestial coordinates of BSSs at |Z| > 17 kpc 
or |Z| < 17 kpc and R > 24 kpc. The BSSs (black circles) and 
BHB stars (red squares) with —200 km s“^ < Vqsr < ~50 km 
s“^ are preferably located over similar Galactic coordinates as 
the Sagittarius Stream. The large dashed rectangles highlight the 
regions where the stream is dominant in (a) The Northern Galactic 
Hemisphere, and in (b) the Southern Galactic Hemisphere. 
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Fig. 19. — Double Residual Root (DRR) plot for the Gaussian 
fits of the distribution of velocities (Vgsr) shown inside the blue 
rectangles of Figure 11, limited by |Z| > 17 kpc. The top panel 
shows the BSSs whose velocities suggest association with the 
Sagittarius Stream (gray shaded region). The bottom panel shows 
the same distribution for the BHB stars. The dashed lines at ± 2 
indicate an approximate 95% significance level. 


shown in Figure 18. This Figure highlights the positions 
(in Galactic coordinates) of BSSs (black circles) and BHB 
stars (red squares) having velocities in the range —200 
km s“^ < Vgsr < —50 km s~^ (the typical Sagittarius 
Stream range of velocities). From inspection of this Fig¬ 
ure, these stars are preferably located over the same coor¬ 
dinates as the Sagittarius Stream (Belokurov et al. 2006; 
Koposov et al. 2012), highlighted by the large dashed 
rectangles: (a) In the Northern Galactic Hemisphere, and 
in (b) the Southern Galactic Hemisphere. The mean he¬ 
liocentric distance distribution found for both the BSSs 
and BHBs with these velocities is Dsg = 21.9 ± 2.2 kpc. 
Within 2(7, these distances agree with the distance to the 
Sagittarius Stream measured by Watkins et al. (2009) us¬ 
ing RR-Lyrm stars: Dsg = 26.1 ± 5.6 kpc. We conclude 
that at least a portion of the held BSSs and BHB stars 
are likely to have originated from the Sagittarius Stream, 
i.e., they are of extragalactic origin. 

Figure 19 shows the excess counts of stars relative to 


the Gaussian hts of the velocity distributions (Vgsr) for 
the BSSs and BHB stars in the blue rectangles of Fig¬ 
ure 11, limited by |Z| > 17 kpc. We have employed a 
variation of a simple difference plot, showing the distribu¬ 
tion of the so-called Double Root Residuals (DRRs; see 
Gebhardt & Beers 1991). The DRR technique takes ad¬ 
vantage of the variance-stabilizing properties of a square- 
root transformation. A DRR plot graphically empha¬ 
sizes where in a distribution there is signihcant lack of 
ht between the data and model - the square-root trans¬ 
formation puts the residuals throughout the ht on an 
equal footing. If the model is an adequate ht to the 
data, then the DRRs are roughly equivalent to normal 
deviates. Thus, a DRR with numerical value exceeding 
±2 is signihcant at the 95% {2a) level, whereas DRRs 
with absolute magnitude less than 1 indicate a reason¬ 
able agreement between data and model. 

The top panel of Figure 19 shows the DRRs for BSSs 
with velocities suggesting association with the Sagittar¬ 
ius Stream (gray shaded region), while the bottom panel 
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shows the same restrictions for BHB stars. There exists 
an excess of about ^400 BSSs and ~250 BHB stars with 
values of Vgsr similar to the Sagittarius Stream. This 
leads to a minimum contamination due to extragalactic 
BSSs arising from the Sagittarius Stream of about 5% 
of the total BSS sample. The contamination for BHB 
stars is even higher, reaching at least 20%. Applying 
the limits described on Figure 11 and the restrictions of 
DRR analysis, we found an excess of 246 BSSs and 146 
BHB stars in the Sagittarius Stream region. Thus, an 
estimate of Fegg/BHE for the Sagittarius Stream is 1.68 
± 0.21, in agreement with the mean value of the rela¬ 
tive blue-straggler frequency in the Galactic halo and in 
the nearby dwarf galaxies within Scr. We also verified 
that the impact of these likely members of the Sagittar¬ 
ius Stream do not change the FBgg/BHB values derived 
by our analysis. 

6. DISCUSSION AND CONCLUSIONS 

We have identified 8001 BSSs and 4796 BHB stars 
from the SDSS/SEGUE spectroscopic database, ob¬ 
tained through DR8. We find that the BSSs become 
more frequent than BHB stars for g > 17, and that 
the mean value for the blue-straggler frequency in the 
Galactic halo is 1.83 ± 0.23, which agrees well with the 
expected upper limit for BSSs in the nearby dwarf galax¬ 
ies (Momany et al. 2007, ~ 2.2) within 2tT. This quan¬ 
tity may be an important input for simulations of stellar 
populations in the Galaxy, in particular for studies that 
include BSSs in these calculations (Chen & Han 2009; 
Conroy & Gunn 2010; Xin et al. 2011; Zhang et al. 2012). 
We also verified that, for the region closest to the Galac¬ 
tic plane that we can study (3 kpc < |Z| < 4 kpc), the 
relative frequency of blue stragglers, Fbss/bhb = 4.11 


± 0.47, agrees well with Preston et al. (1994) BSSs fre¬ 
quency estimate in the Solar Neighborhood. The relative 
frequency drops with distance from the plane, reaching 
^ 1.5 — 2.0 in the inner-halo region; this ratio continues 
to decline to ~ 1.0 in the outer-halo region. 

Another interesting result from our analysis is the ver¬ 
ification of the extragalactic origin of a small fraction 
(^ 5%) of BSSs (a somewhat higher fraction, ^ 20%, is 
obtained for BHB stars) found among stars in the outer- 
halo region. 

The distribution of BSSs throughout the Galaxy can be 
used to constrain how the frequency of primordial binary 
systems (assuming this is the primary formation channel) 
varies with distance from the Galactic center. Our anal¬ 
ysis indicates that the relative number of BSSs decreases 
with both R and | Z |, thus it would follow that the relative 
frequency of primordial binary systems should decrease 
as well. 
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